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The dc Josephson effect refers to the dissipationless electrical current - the supercurrent 
- that can be sustained across a weak link connecting two bulk superconductors.^ This ef- 
fect is a probe of the fundamental nature of the superconducting state, which depends cru- 
cially on the spatial dimensionality of the superconducting electrodes. For bulk (i.e. three- 
dimensional) superconductors, the superconductivity is most robust, and the Josephson ef- 
fect is sustained even at nonzero temperature. However, in wires and thin- films {i.e. lower- 
dimensional superconductors), thermal and quantum fluctuations play a crucial role. In 
superconducting wires, these effects qualitatively modify the electrical transport across a 
weak link.^ Despite a number of recent experiments involving weak links between thin-film 
superconductors,^"^ little theoretical attention has been paid to the nature of the electrical 
conduction in such systems. Here, we analyze the case of two superconducting thin films 
connected by a point contact. Remarkably, the Josephson effect is absent at nonzero tem- 
perature, and the resistance across the contact is nonzero. Moreover, the point contact re- 
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sistance is found to vary with temperature in a nearly activated fashion, with a universal 
energy barrier determined only by the superfluid stiffness characterizing the films, an angle 
characterizing the geometry, and whether or not the Coulomb interaction between Cooper 
pairs is screened. This behavior reflects the subtle nature of the superconductivity in two- 
dimensional thin films, and should be testable in detail by future experiments. 

Soon after the development of the microscopic theory of superconductivity,' Josephson pre- 
dicted a remarkable manifestation of it: charge is transferred via a supercurrent across a weak 
link (e.g. an insulating barrier) between two bulk superconductors, even in the absence of a volt- 
age difference.^ This remarkable prediction by Josephson was already observed at the time by 
Giaever,^ but was misinterpreted as a "metallic short."^ Both were awarded the Nobel prize in 
1973 for their discoveries. 

Since its discovery, the Josephson effect has had a major impact on a broad spectrum of tech- 
nologies. The most sensitive magnetic flux and electromagnetic radiation detectors are SQUIDs 
(Superconducting Quantum Interference Devices), consisting of two Josephson junctions con- 
nected in parallel. SQUIDs play a crucial role in many condensed matter experiments and in radio- 
astronomy, and recently in biomagnetic detectors to monitor brain activity. Josephson junctions 
were also a stepping stone in computer technology. In the 1970's, the burgeoning semiconductor 
revolution led IBM to abandon its effort to construct computers out of arrays of Josephson junc- 
tions. But today, circuits made of Josephson junctions are back in vogue, in the effort to construct 
a quantum computer. ^"^'^^ 
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Josephson's prediction was based on the fundamental principle of broken symmetry. A su- 
perconductor can be thought of in terms of a complex- valued wave-function (or "order parameter") 
?/^(r) for Cooper pairs. At the onset of superconductivity, Cooper pairs undergo Bose condensation, 
and the order parameter becomes nonzero with a well-defined phase Lp: ipir) = \^\ exp{i(p). The 
symmetry of phase rotations is thus broken. When the phase (/?(r) varies in space, a bulk super- 
current proportional to its gradient results. Similarly, a Josephson supercurrent flows between two 
superconductors with different phases. It is proportional to the sine of the phase difference across 
the weak link: 

^ = sin ((/?!- (y^s) . (1) 

n 

Here, J is the Josephson coupling energy, which characterizes the coupling strength between the 
superconductors. 

The Josephson effect was initially measured between two bulk (i.e. three-dimensional) su- 
perconducting electrodes. Experimental technique has since evolved dramatically. Present-day 
experiments probe superconductivity in nanometer-scale samples, as well as in systems with re- 
duced dimensionality. Particularly intriguing recent experiments measured the electrical resistance 
across nanowire junctions between thin-film superconductors,^ '* through a narrow constriction 
(about 20 nm wide) between two films as thin as 2.5 nm,^ and between a film and a superconduct- 
ing scanning tunneling tip.^ In such nanoscale systems with reduced dimensionality, the simple 
Josephson effect of Eq. ([T]) must be revisited. 

In large bulk superconductors, the phase of the order parameter is a rigid, essentially 
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classical, variable. In equilibrium, and in the absence of currents, ip is locked to a single fixed 
value; the system exhibits long-range order (LRO). This is true even in the presence of thermal 
fluctuations, provided the sample is below the superconducting transition temperature. The phase- 
rigidity of three-dimensional (3d) superconductors prevents phase fluctuations that would suppress 
the Josephson supercurrent, and Eq. ([T]) holds. It is important to emphasize that these statements 
hold in the limit of large superconducting electrodes; otherwise, there is neither true LRO nor a 
true Josephson effect. 

In thin-film superconductors, thermal fluctuations of the phase dramatically alter this simple 
picture. Specifically, phase correlations are no longer infinitely long-ranged; instead, the corre- 
lator ^e"^*^^^^e~*'^*^^2)) (where (...) denotes a thermal average) decays as a power law with spatial 
separation. This phenomenon is known as quasi-long-range order (QLRO), and occurs below 
the Berezinskii-Kosterlitz-Thouless phase transition temperature Tbkt- Thermal phase fluctua- 
tions varying smoothly in space, present due to thermal excitation of the superconducting plasmon 
mode, are responsible for the power-law phase correlations. Only at T = 0, where these are 
frozen out, does LRO obtain. Remarkably the resistance of a two-dimensional (2d) film vanishes 
for T < Tbkt despite the absence of LRO. Above Tbkt, however, the superconducting state is 
disrupted by topological vortex defects, around which Lp{r) winds by an integer multiple of 2tx. A 
nonzero density of mobile vortices scrambles the phase; only short-ranged superconducting corre- 
lations survive, and the above correlator decays exponentially with spatial separation. In addition, 
the resistance of the film is nonzero in this regime. Below Tbkt, isolated vortices are expelled 
from the film or tightly bound into vortex-antivortex pairs, and QLRO results. Below, we consider 
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the fate of the Josephson effect across a point contact separating two films below Tbkt- 

The question of the Josephson effect in systems with reduced dimensionality has been ad- 
dressed previously, but in the context of one-dimensional (Id) superconducting wires. At any 
nonzero temperature, the phase correlations in Id are short-ranged, and one expects a non- vanishing 
resistance even in the absence of a weak link.^^'^^ This can be understood in terms of thermally 
excited phase-slip events. At zero temperature, a Id superconductor exhibits QLRO, and the re- 
sistivity vanishes. This holds provided g > Qc, where g is the superfluid stiffness measured in 
appropriate units, and is proportional to the cross-section of the wire. Qc is on the order of unity 
and is determined by the nature of any Umklapp scattering that may be present. The effects of a 
weak link on a Id superconductor have been addressed theoretically; provided g > was found 
that at zero temperature the weak link "heals" itself, and the resistance across the point contact 
vanishes (except in the case of a wire with only a single transverse electron mode at the Fermi en- 
ergy). At small nonzero temperatures, the resistance across the point contact is predicted to vanish 
as a power law in temperature:^ 

R(xT^^s-^\ (2) 

The Josephson effect is obliterated at nonzero temperature, in dramatic contrast to the weak link 
between 3d superconductors. (It should be noted that, although ref. 2 explicitly dealt with a single- 
channel quantum wire, the effective field theory is identical to that for a many-channel supercon- 
ducting wire, and the preceding statements follow immediately from the analysis there.) 

What is the fate of the Josephson effect in the intermediate case, when the point contact 
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is between two films rather than wires or bulk electrodes? Surprisingly, despite several recent 
experiments on 2d films, very little theoretical attention has been paid to this problem (but see 
ref. 14). In this Letter we determine the tunneling resistance across a point contact between two 
films in the geometry shown in Fig. 1. As it turns out, a weak link between two films is almost 
superconducting, but a true Josephson effect is absent except at zero temperature. Specifically, we 
find that the resistance R{T) drops very rapidly upon cooling, in a nearly activated fashion. At low 
temperatures 



R{T) tl 1 r EAiT) 



exp 



(3) 



where ty is an amplitude for quantum phase-slip processes, discussed in more detail below. This 
formula is expected to be asymptotically exact at very low temperatures (see below). Here Rq — 
h/Ae^ is the quantum of resistance, and Ea{T) is a temperature-dependent activation energy: 

Ea{T) = + {2cKjTr^jy ^"^^ 

R{T) is positive for all nonzero temperatures, but vanishes in an activated fashion as T — > 0, up 
to the logarithmic correction contained in the denominator of Ea{T). R{T) thus vanishes faster 
than any power law in the low-temperature limit. The scale of the activation energy is set by the 
superfluid stiffness in the 2d films. Kg — h^ng/m, where Ug is the density of Cooper pairs in the 
film and m is the pair mass. The dimensionless number c is given by 

where 9 is the "opening angle" shown in Fig. 1, and a is a parameter characterizing the range 
of the interactions. Specifically, a = 2 for Coulomb interactions and a — 1 for screened in- 
teractions, as can be obtained in the presence of a superconducting ground plane (see below). 



Remarkably, the main effect of the Josephson coupling is to set a crossover temperature Tj ~ 
{hujc/2kB)ex.\){—2cKs/'K'^J), where uj^, is a high energy cut-off discussed further below. At low 
temperatures, where the ln(l/T) dominates the denominator of Ea{T), i.e. T <ti Tj, the depen- 
dence on J disappears altogether, and one recovers the universal result 

Ea ~ cKj\n{nwj2kBT). (6) 

The important temperature scale in this limit is set hy Kg, a property of the 2d films that can be 
measured independently^"^ and is independent of the details of the contact. At higher temperatures, 
R{T) is approximately given by a purely activated form {i.e. no logarithmic corrections), with the 
barrier height set by the Josephson coupling: Ea ~ ^J- 

The resistance formula of Eq. ^ is exact at low temperature in the sense that 
ln[_R(r)]/ ln(_Rineasured) — 1 as T — > 0. The ratio of the resistances themselves does not go to 
unity, as there are additive corrections to Ea{T) proportional to ln~^(l/T). These corrections are 
contained in the integrals leading to Eq. © (supplementary information), so one can do better by 
evaluating them exactly, which must be done numerically. 

Our result for R{T) fulfills some basic physical requirements. The zero-temperature resis- 
tance vanishes, as it should in the presence of the true LRO that obtains in the 2d films. Further- 
more, the resistance vanishes faster than any power law as T approaches zero, as is reasonable 
upon comparison with the Id and 3d cases. Finally, the activation energy increases monotonically 
with increasing Josephson coupling. 
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Before addressing the experimental implications of our result, we briefly discuss its deriva- 
tion and the underlying physics. The thin-film superconducting electrodes and the point contact are 
modeled by a quantum phase Hamiltonian, which focuses on the quantum and thermal fluctuations 
of (p{r). This is legitimate at temperatures well below the quasiparticle gap scale and Tbkt- The 
Hamiltonian for a single thin-film electrode is, 



where n{r) is the fluctuating Cooper-pair density, canonically-conjugate to (p{r). The spatial inte- 
grations range over the area of the film. The first term encodes the energy cost for phase gradients; 
this is essentially the kinetic energy of the superflow. The second term is a density-density inter- 
action; we consider both Coulomb [V{r) = (2e)^/|r|] and screened [V{r) = {fi?v1/ Ks)5{r — r')] 
interactions. This Hamiltonian describes the quantum dynamics of the plasmon mode of the su- 
perconductor. In the case of screened interactions, there is a linearly dispersing acoustic plasmon 
(a; = Vgk), whereas for Coulomb interactions the plasmon has the dispersion huj = y/Sire'^Ksk. 
The high-energy cutoff lUc appearing in Eq. Q is given in terms of this dispersion: cUc = uj{k = 
2n/^), where ,^ is a short-distance cutoff on the order of the superconducting coherence length. 

The second element in our system is the weak link, which we shall refer to as a point contact. 
We model it as a point-like Josephson coupling between the two electrodes, and set the coordinate 
origin for both films (i.e. r = 0) at the contact: 



Here fi,2{r) are the phase fields in the two electrodes. This term can be interpreted as a process 
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contact 



Jcos(y.2(0)-</^i(0)). 



(8) 
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hopping Cooper pairs across the contact, where J is the hopping ampUtude. It should be noted 
that we have neglected inter-film Coulomb interactions. Because any geometric, zero-frequency 
capacitance will diverge with the size of the films. Coulomb blockade effects will be unimportant 
except below an extremely low temperature. Above this temperature, we do not expect inter-film 
interactions to modify the low-temperature resistance. 

In the limit J <^ Kg (i.e. a poor contact), we may calculate the current response to a small 
voltage bias in terms of the hopping of Cooper pairs across the contact. Specifically, one can 
expand in powers of J to calculate the conductance G across the point contact: G{T) = 02 + 
0( J^). Kim and Wen^^ found that the coefficient a2{T) is divergent for temperatures below T* — 
9Ks/2akB, and concluded that the conductance diverges (and hence the resistance vanishes) for 
T < T*. However, the physical meaning of this result is unclear; in fact, we shall show that it 
indicates not zero resistance but rather a breakdown of perturbation theory. More information, and 
an independent test of this perturbative result, can be obtained via an expansion in the opposite 
limit: that of a good contact. Indeed, the above result already suggests a strong tendency towards 
this limit, and justifies the approach taken below. 

The limit of a good contact is readily accessed through a dual picture. Here, we begin 
by assuming perfect phase-coherence across the contact. This coherence is then weakened by 
quantum phase-slip events, which occur when a vortex tunnels across the film near the contact, in 
the direction transverse to the current flow (as illustrated in Fig. 1). A phase slip makes the relative 
phase of the electrodes wind by 2%; this 2% twist can either heal locally by a phase slip in the 
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opposite direction, or propagate outward into the electrodes. The latter process will register as a 
voltage spike across the system, due to the Josephson relation 

In this dual picture the resistance can be obtained through a perturbation expansion in t^, 
which is the amplitude for a phase slip to occur, or, equivalently, for a vortex to hop through the 
contact. Formally, we integrate out the degrees of freedom in the films to obtain an action for the 
phase difference = <fi{0) — ^2(0) across the contact. The dual action for the phase slips can 
then be obtained by a Villain transformation.^^ The phase slips are instanton events where jumps 
by ±27r, and we define as their fugacity. The single vortex hopping process with amplitude ty 
should be viewed as an encapsulation of the many different physical vortex-hopping processes in 
the vicinity of the contact (see Fig. [T]). We note that ty depends only weakly on temperature (it 
goes to a constant at T = 0). It also has an implicit dependence on J, decreasing as J increases. It 
would be interesting to compute this dependence, but we do not attempt to do so here. 

Following the above approach, we have calculated the voltage response to a vanishingly 
small current across the point contact, obtaining the resistance formula Eq. Q at order tl- This 
formula encapsulates the exact low-temperature dependence of the tl term, and is quantitatively 
accurate over a reasonably large temperature range. As an example, using the parameters for the 
system of MoGe films with ground plane discussed below, one can compare -R(T) to the exact 
resistance at order tl (obtained numerically). The logarithms of these quantities agree within 15% 
for T < 4K. Furthermore, contributions of higher order in ty involve larger effective activation 

10 



barriers and are unimportant at low temperatures. This is to be expected, as these terms correspond 
to processes involving more phase-slip events, and the phase slips are strongly suppressed by the 
QLRO in the films. 

The resistance formula obtained from the dual calculation is expected to be exact at low 
temperatures, even in the case of a poor contact. Physically, the expansion in phase-slip events is 
valid in the limit of small resistance. As R{T = 0) vanishes due to the presence of LRO in the 
films, the phase-slip expansion is always correct at low temperatures. Reasoning along the same 
lines, our result should be a good approximation as long as the argument of the exponential in 
Eq. 01) is large and negative. 

The most interesting regime is the universal low-temperature limit, i.e. Eq. This regime 
will be more easily accessible when J/ Kg is rather large; this could be achieved by fabricating the 
point contact as a short and wide constriction. 

The most important quantity determining the universal low-temperature resistance is the su- 
perfluid stiffness Kg. This can be obtained directly for thin films by an inductance measurement. 
In fact, one could make a direct check of our theory by varying Kg in situ with an in-plane mag- 
netic field. The other relevant parameters are the Josephson coupling J and the cutoff frequency 
ujc, which should be considered fitting parameters. In many cases it should be possible to estimate 
J from Kg and the geometry of the contact (see below). An estimate of ujc is given in terms of the 
coherence length in the films, as discussed above. 
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Our result should be accessible to experimental tests in a variety of different systems. Here 
we consider the specific case, which was already realized by Chu et al, of a narrow constriction 
between two MoGe films.^ We imagine a modification of their setup, in which a superconduct- 
ing ground plane with stiffness is added a distance d below the 2d films, which we take to be 
100 nm. We consider ^ K^, as appropriate for a thick ground plane, to screen out the Coulomb 
interactions as much as possible. This setup leads to an acoustic plasmon propagating in the films 
with velocity Vg = ^^/IGne'^Ksd. There is also a plasmon with u cx \/k dispersion, which propa- 
gates within the ground plane and does not contribute to the tunneling transport. The constrictions 
in the experiments of ref. 6 had width w ~ 20 nm and length ^ 100 nm. The films themselves 
exhibit a superconducting transition at Tbkt ~ 5 K, and have a coherence length ^ ~ 7 nm. We 
consider 9 — n and would have a — 1 with the screening ground plane in place. For illustrative 
purposes, we assume Ks/ks ~ lOK for T <C Ibkt, and then huc/kB ~ 27rhvs/kBd 2000K. 
(In the presence of the ground plane, the short-distance cutoff in the films is set by d rather than 
^.) We estimate J/kB ^ Kgw/ikB ~ 2K. For these parameters, neglecting the temperature- 
dependence of Kg, the argument of the exponential in R{T) is —1 for T pa 5K, so we conclude 
that our formula should be good for T < 5K. 

One of the most important considerations for experiments is that the superconducting films 
be in the thermodynamic limit. For a film of linear dimension L, we can define a crossover tem- 
perature scale /cbTl = hu;{k — 27r/L), which vanishes as the linear size L of the films diverges. 
When T ^ Tl, thermal processes prevent quanta of the plasmon mode from traveling coherently 
between the film edges. This thermal decoherence washes out the influence of the environment 
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(and other finite-size effects) on the system of films and point contact. With a superconducting 
ground plane as above, one has ksT^ = 2'Khvs/ L. For L = 1cm we find Ti ^ 2 x 10~^ K. This 
leaves a broad temperature range over which the resistance formula Eq. Q should apply, and also 
points out that it is important to have rather large films to lower T^. It is difficult to achieve such a 
broad temperature range without the ground plane, as Ti is only proportional to L^^/^ in that case; 
for the parameters as above, but with no ground plane, ^ IK. However, surrounding the film 
by a medium with rather large dielectric constant e would lead to a modest improvement, because 
Tl oc 

It is interesting to point out that the system studied here cannot be understood in terms of the 
phenomenological RCSJ (resistively and capacitively shunted junction) model, which is one of the 
standard theoretical approaches to Josephson junction physics."' In that model, one imagines that 
the junction is shunted by a parallel resistor and capacitor. One then writes down a Hamiltonian in 
terms of only 0, the phase difference across the junction, and N, a canonically conjugate integer- 
valued Cooper pair density. This takes the form: 

Hrcsj = - Jcos(0), (10) 

where C is the capacitance. The resistor is a phenomenological model for the dissipation in the 
system, which is assumed to be of the Caldeira-Leggett (CL) form'^ appropriate for dissipation due 
to the plasmons of a Id superconductor,^ or the similar Ambegaokar-Eckem-Schon (AES) form^^ 
that models dissipation due to a Fermi surface of gapless quasiparticles or electrons. The RCSJ 
model does not apply, however, for point contacts between large 2d or 3d superconducting leads 
with gapped quasiparticles. In these cases, the capacitance diverges as the size of the leads is taken 
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to infinity. The charging energy vanishes in this Umit, and Coulomb blockade effects disappear 
above a temperature scale proportional to 1/L for both 2d and 3d cases with Coulomb interactions. 
As for the resistor, in these systems there is either dissipation due to the superconducting plasmons 
in 2d, or a pinning effect due to LRO in 3d (supplementary information). These effects cannot 
be modeled by dissipation of the CL or AES form, as they are much stronger at low frequency. 
Roughly speaking, the effect of the 2d or 3d electrodes is much closer to connecting the junction to 
a pair of ground planes, rather than to a resistor. Indeed, even if one adds a physical shunt resistor 
to the system, at low frequencies it will be shorted out by the ground planes. 

Finally, we expect that our result will describe transport in a variety of systems - our model 
only requires that the low-energy excitations in the films can be described by some bosonic degrees 
of freedom governed by an XY model in its QLRO phase. This physics is most readily accessible 
in 2d superconducting films, and we hope our work will stimulate more experiments on these 
intriguing systems. 
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Supplementary Information Here, we provide further information on the technical aspects of our results. 
We shall discuss the action for the phase difference cj) across the point contact, the duality transformation to 
the phase-slip action, and the derivation of the resistance formula. In this supplementary information, we 
work in units where h = \ and ks = 1- 

The partition function for the coupled system of films and point contact can be written as the phase-space 
functional integral 

Z = J 'D(fi{r,T)T>ni{r,T)ex.p {- Si - S2- Sj), (11) 

where 

Si=i dr d\dripi{r,r)ni{r,T)+ / (irH^i^(r) (12) 
Jo J Jo 

is the action for film i (i = 1,2), and 

Sj = -J drcos(v72(0,r)-v?i(0,r)) (13) 
Jo 

is the action for the point-contact Josephson coupling. A product over space, time and film indices is implicit 
in the integration measure; that is, 

V^i{r,T) = J] \{d^i{r,T), (14) 

j=l,2 r,T 

and similarly for 'Dni{r, r). The Hamiltonian Wgjj^ is defined in Eq. ©, and (3 = 1/T. The fields ipi{r, r) 
and ni{r,T) satisfy periodic boundary conditions in imaginary time, e.g. (pi{r,T + /?) = ipi{r,T). Fur- 
thermore, at the edges of each film, the phase satisfies the boundary condition V(/?j • = 0, where v is the 
normal vector to the edge. This is simply the physical requirement that no supercurrent flow across the edge. 

Because the films are governed by a quadratic action, standard techniques can be employed to integrate out 
all the fields except (I){t) = (/?2(0, t) — (/5i(0, r). In order to do this, it is first convenient to transform the 
fields to a basis in which Si is diagonal. In the case of short-ranged interactions {i.e. V{r — r') oc 5{r — r')) 
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this is entirely straightforward; the kinetic energy term is diagonal in a basis of eigenfunctions of the Laplace 
operator — with the above boundary condition, and in this case the density-density interaction is diagonal 
in any basis. In the case of an infinite, translation-invariant plane [i.e. 9 = 27r], Coulomb interactions can 
also be treated in an entirely straightforward manner. In general. Coulomb interactions are somewhat more 
difficult to deal with, because it is nontrivial to simultaneously diagonahze the kinetic and interaction terms. 
To avoid this difficulty, we replace the Coulomb interaction with an effective potential, which is consttucted 
to satisfy certain key properties that we beheve are enough to leave the final result unaffected.^ 

At this point, a straightforward calculation yields an effective partition function depending only on 0: 

Zes = J V<I){t) exp {-So- Sj), (15) 



where 



l</'K)r. (16) 



Here, ujc is a high-frequency cutoff set by the short-distance cutoff in the films, ujn = 2iml(3 (n integer) are 
the usual Matsubara frequencies, and 

0(a;„)= f dTe'^-^(t>{T). (17) 

JO 



^Consider an infinite, translation-invariant plane. There, the Coulomb interaction is diagonal in momentum space, 
and has the Fourier transform V{k) = Sttc^ /\k\. The Laplacian is diagonal in a plane- wave basis, and we can define 
its square root by writing V— V^e'*^'' = \k\e^^'^. Therefore we can write = 87re^(v'— V^)~^. The new potential 
V" (r, r' ) is constructed so that this equation also holds, except that V— is defined in terms of eigenfunctions of — 
satisfying the V(/? • w = boundary condition. The resulting form depends separately on r and r', and not only on 
the distance |r — r'| ; however, as the problem in the wedge geometry already lacks translation and rotation invariance, 
it should be harmless to sacrifice these in the form of the potential, as far as universal properties are concerned. 
Furthermore, and most importantly, it can be shown that V has the same scaling behavior as the Coulomb potential: 
V{sr, sr') = \V{r, r'), where s is a positive real number. 
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It is interesting to compare this result to its 3d analog. There one begins with the Hamiltonian for the bulk 
electrode 



_ -^s 



J (fr{Vipf + d^rd^r'n{r)n{r')V{r - r'), (18) 



2 

where Kg"^ is the 3d superfluid stiffness in units of energy per unit length. Proceeding as above, one finds 

iy-3d 

where g{u!n) goes to a nonuniversal constant as a;„ 0. This is a "mass term" for 0, which localizes 
the phase and is expected to lead to a tunneling resistance of strictly zero (for superconducting leads in the 
thermodynamic Umit). This expectation is confirmed by a calculation analogous to that described below for 
the 2d case. 

Furthermore, it is instructive to compare Sq and Sq'^ to the Caldeira-Leggett term that would be included in 
an RCSJ model. This takes the form 

5o°^cx^K||.^(a;„)|^ (20) 

which is clearly dominated by both and Sq"^ at low frequency. Therefore, even if it were appropriate to 
include a term of the form Sq^ (it is not for the systems discussed here), it would have no effect on the 
low-temperature physics. 

Now we return to the 2d system. The partition function Z^s is the starting point for the duaUty transformation 
that allows us to work directly in terms of phase sUps. The first step is to replace the cosine with a Villain 
function: 

oo J 

exp(Jcos(0)) ^ ^ exp (--((/.- 27rr?)2) (21) 

7J= — oo 

This is simply a different 27r-periodic potential, the detailed form of which is not expected to affect universal 
properties. Strictly speaking, when making this replacement one should also replace J by an effective 
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Josephson coupling (which is expected to be approximately equal to J). However, because J only enters 
the final result as a fitting parameter, it is not necessary to make this distinction here. 

We obtain the Villain partition function by making r] a periodic function of imaginary time (ry(r + /3) = 
r/(r)): 

r ^ r J r0 . 

(22) 



Zv = j V(t) ^ exp 

»7(t)=— oo 



^0-^^ dT{<l^{T)-2TTrj{T)f 



This can be specified more precisely by writing r/(T) = dr' p{t') and working in terms of p{t), which 
can be interpreted as the imaginary-time density of phase slips. (We can set t/(0) = without loss of 
generality.) Due to the periodic boundary conditions, there are an equal number of "positive" {rj ^ rj + 1) 
and "negative" (r/ — > r/ — 1) phase slips in each configuration summed over to form Zy- Letting m be the 
number of positive phase slips in a given configuration, we can write 



P(r) = X;fe-r,)-<5(T-f,)), 



(23) 



where Tj and are the times at which positive and negative phase slips occur, respectively. In Zy, the sum 
over T] becomes a sum over p{t), which can be defined precisely as 



f 



j dTi... dTm j dri. . . dT„ 



(24) 



p(t) m=0 

Note that the phase-slip fugacity has now explicitly appeared. Upon transforming to Fourier space and 
integrating out (f), we obtain a partition function only in terms of p: 



Zv = ^ exp 

P{r) 



-y 

2/3 ^ 



J 



2T^p{iO„ 



(25) 



The duality transformation is completed by observing that Eq. d25l is equivalent to the following sine- 
Gordon theory: 

^duai = j Ve{T) exp ( - 5;?""' + U dT cos(27r0)) , (26) 
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where 



cdual 



1 r 1 vr 



J AcKs 



In 



(27) 



Expanding in powers of ty we obtain 



n=0 ■ a-i---(T„=±l i=l 



(28) 



At this point it is helpful to note that 9{t) — > {0{t) + constant) is a symmetry of S'q"^', which means 
that only "neutral" configurations with J2i = will contribute. This allows us to reorganize the above 
expression and write: 



^duai = Y,Jveexp[- S^'^'^'ie] + 2m J dTp{T)9{T) 



(29) 



Pir) 

Upon integrating out 6 we immediately obtain Eq. d25t . 

These manipulations also show that the operator exp(27ri0(r)) inserts a phase slip at time r. As a result, the 
voltage across the contact (in a Hamiltonian formulation) is 



y = ^sin(27r0). 
2^ 



(30) 



The linear response resistance (at leading order in t^,) can be expressed in terms of a dual Kubo formula: 



R(T) = - lim Re 



-G)j(^) 

LO 



(31) 



where 



r 

Jo 



die— *([y(t),y(0)])i„=o. 



(32) 



Note that the thermal average is evaluated in the limit ty = 0, i.e. in the quadratic theory governed by Sq^^^. 
This Green function can be obtained by analytic continuation from imaginary time, resulting in an integral 
expression for R{T). It is convenient to deform the contour in this integral, leading to the following simple 
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form: 



R{T) = ^I dtef^'\ (33) 



where 

dx cos(xt) — cosh(^) 



fit) = IGcKs / ^ ^ — o. (34) 



The function /(t) has a maximum at t = 0, and /(O) ^ — oo as T ^ 0. Therefore, the saddle-point method 
can be used to evaluate R{T) at low temperature; one finds 



This formula is our main result, and its evaluation at low temperature leads to Eq. In particular, at 
asymptotically low temperature 

rK 

The full form of Ea{T) given in Eq. Q is an approximation to /(O) that reduces to the universal form at 
low temperature, and is valid over a larger range of temperatures. For example, in the system of MoGe films 
in the presence of a superconducting ground plane discussed in the paper, Ea{T) is within 15% of its exact 
value [obtained by a straightforward numerical evaluation of /(O)] for T < 4K. On the other hand, the 
universal limit Eq. (I36t only agrees within 50% in this temperature range. 
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Figure 1: Illustration of point-contact Josephson tunneling between two thin-film superconduc- 
tors (shaded regions), in the wedge geometry characterized by the angle 0. This geometry directly 
models the experimental setup of ref. 6. Charge transport across the contact can be viewed most di- 
rectly in terms of the hopping of Cooper pairs (solid arrow) with amplitude given by the Josephson 
coupling energy J . Alternatively, one can consider phase-slip events, in which the phase difference 
between the two films winds by 27r. These events occur when a vortex hops across one of the films 
near the junction, transverse to the current flow (dash-dot line); this produces a momentary voltage 
spike between the films. Because vortices can be ignored within the films at low temperature, all 
such events can be lumped together into a single effective process where a vortex tunnels directly 
across the contact (dashed line) with amplitude t^. As we argue in the text, it is essential to adopt 
the phase-slip point of view to understand the low-temperature transport. 
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Spatial dimensionality 


Resistivity of individual electrodes 


Point contact tunneling resistance 


d = 3 


Zero 


Zero 




Zero 


Nearly-activated 


d= 1 


Power- Law 


Power-Law 



Table 1 : Comparison of the electrode resistivity, and tunneling resistance across a point 
contact, for superconductors of varying spatial dimensionality. The entries of the table de- 
scribe the behavior at low temperature. The main result of our work, which completes this 
table, is the nearly-activated behavior of the resistance across a point contact between 
two 2d superconductors. 
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